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ABSTRACT 

This paper discusses the effects of variations in 
system parameters on the dynamic behavior of a Free- 
Piston Stirling Engine/Linear Alternator (FPSE/LA)-load 
system. The mathematical formulations incorporate both 
the mechanical and thermodynamic properties of the 
FPSE, as well as the electrical equations of the con- 
nected load. State-space technique in the frequency 
domain Is applied to the resulting system of equations to 
facilitate the evaluation of parametric impacts on the 
system dynamic stability. Also Included is a discussion 
on the system transient stability as affected by sudden 
changes In some key operating conditions. 

Some representative results are correlated with 
experimental data to verify the model and analytic for- 
mulation accuracies. Guidelines are given for ranges of 
the system parameters which will ensure an overall stable 
operation. 

INTRODUCTION 

The free-piston Stirling engine/linear alternator 
(FPSE/LA) is an attractive, thermoelectric energy conver- 
sion system for space applications due to its potential for 
long life and reliability. However, when integrated into a 
power system with connected load(s), system perfor- 
mance must be evaluated to ensure that unfavorable 
system interactions do not occur. Treatment of such 
interactions is not commonly found in the existing litera- 
ture which deals mainly with establishing the engine 
operating conditions (Refs.1 to 7). Additionally, a combi- 
nation of detailed representations of the eng ine thermody- 
namics and nonlinearities, and load modeling has not 
received widespread attention. 

This paper describes a general dynamic analysis of 
a FPSE/LA connected to a load. The objective of the 
paper is two-fold. First, the paper demonstrates the 
application of the state-space technique (Ref. 8) to model 
the dynamics of a FPSE/LA-load system. This technique 


is then used to predict the impact of some key system 
parameters and operating conditions on the dynamic 
stability of a Space Power Research Engine (SPRE)/LA- 
load system at NASA Lewis Research Center (Ref. 9). 

The mathematical model used for the FPSE Is 
based on the recommendations in Ref. 10. The analysis 
incorporates the engine dynamics and thermodynamics, 
and detailed modeling of the LA and connected load(s). 
The state-space formulation used in the analysis is a 
modern control theory approach of representing a system 
by its state variables. By definition, a minimum set of n 
state variables, typically denoted by *X", Is necessary for 
a complete description of the internal status or state of 
the system. The state-space designates the n-variable 
co-ordinate space in which X ranges. The representation 
of the system equations by the compact state-space 
model (SSM) permits the application of powerful vector- 
matrix theory, and can readily be relegated to a computer 
to yield a complete description of the system. 

The simulation, performed on the MATLAB software 
(Ref. 11), yields results in the form of root locus plots of 
the system eigenvalues (Ref. 8). The plots show the 
migration of the system roots caused by parametric 
variations and, hence, the effects of such variations on 
the stability margins of the SPRE/LA-ioad system. The 
advantage of the frequency-domain-based eigenvalue 
analysis is Its single computation of the exact modes of 
oscillation of the dynamic system. Also, it can comple- 
ment, valuably, Information obtained by time-domain 
simulation and testing of a physical system. Many 
practical systems are nonlinear. However, the usefulness 
of the linearized eigenvalue analysis Is Its use In identify- 
ing potential worst-case operating conditions which can 
then be verified by nonlinear time-domain analysis and, 
possibly, testing. Hence, by using the appropriate input 
data, the analytic method and results such as reported in 
this paper can be used to expedite, or assist in, the 
development of requirements for the application of a 
FPSE/LA In a space power system. 


STUDY SYSTEM 

The analytic formulations discussed in this paper 
relate to a system which consists of a single cylinder 
FPSE connected to an LA via the engine power piston. 
Figure 1 depicts a block diagram representation of the 
study system. The LA output terminals in series with a 
tuning capacitor, C T , feed through a load controller to an 
external load. The controller may be used to "dump" 
excess generated electric power, in the case where 
power demand (by the external load) is less than the 
power supplied by the engine/alternator source. In this 
case, the load controller acts to maintain power balance. 
When the load demand exceeds the generated power, 
the engine will tend to stall. 

ANALYTIC APPROACH 

The overall performance equations of the FPSE/LA- 
load system are developed in three stages. First, the 
equations describing the dynamics of the FPSE are de- 
rived. This is followed by the subsystem comprising the 
LA, parasitic load and an external load. Finally, the 
combined system equation is formulated for subsequent 
analysis. 


The element C DD is the total damping coefficient of 
the displacer; that is, the sum of the displacer -to -engine 
casing (or ground), C 0 , and the displacer-piston coupling 
C PD . Similarly, the total piston damping C pp comprises 
the self term C p and piston-displacer coupling C pp . 
These definitions also apply to the stiffness matrix, |KJ. 
The external force vector [F E ] consists of the displacer 
and piston components F E0 and F EP , respectively 
Typically, F eo is nonexistent. 

The matrices |C T J and [K T ] are associated with 
the engine thermodynamic force vector [F T ] expressed 
in EQ 2. The negative sign denotes the restoring nature 
of the force [F T ]. 

I'M * -l|C r ]|X] - iK,m (2) 

The details of derivations of the engine equations are 
documented in Ref. 12. The force [F T ] is expressed in 
EQ (3), where A„, A„ and A„ are the displacer, displa- 
cer rod and piston areas, respectively. The parameter P c 
is the compression space pressure on the displacer and 
piston. The term P HX is the heat exchanger pressure 
difference between the expansion pressure P E and P c . 
The displacer and piston components of [F T ] are F TD 
and F^,, respectively. 


DYNAMIC EQUATIONS OF FPSE 

The following simplifying assumptions are made - 

(1) Schmidt’s thermodynamic analysis Is evoked such 
that: 

(a) the displacer and piston motions and the 
working space pressure are sinusoidal 

(b) the working fluid obeys the ideal gas law, and 
expands and compresses isothermally 

(2) The working space gas pressure is constant but time- 
variant. 

(3) The bounce space pressure balances the average 
working space pressure. Hence, average positions of the 
power piston, displacer and cylinder casing are stationary. 

The dynamic equations of the FPSE/LA system, In- 
corporating the mechanical and thermodynamic properties 
of the Stirling engine, Is stated In EQ (1). 

[M m ][XJ * [C M * C T 1(X] * (K„ • K t 1[X| . (F t ) (D 


where the mechanical system matrices are: 


[MJ • displacer (D) and pistion (P) mass matrix * 
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The thermodynamic pressure P T is assumed a linearly 
dependent function mainly of position and velocity. 
Hence, EQ (3) may be recast into EQ (4), where the 
superscript P denotes partial derivatives of the pressure 
terms in EQ (5). 
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and m D and m p are the displacer and piston masses, 
respectively; 


and 


[CJ « damping matrix - 






Comparison of EQS (2) and (4) implies the relationship of 
EQ (5). 


icj - -IAiJICT < 6 > 
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JKJ - stiffness matrix 
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Equation (5) requires that each of and P c must be 
a continuous function of the piston position and velocity. 


EVALUATION OF [K P| AND [C p ) 


The matrices [K^ and [Cj and, hence, [Ky] and [C T ] 
may be obtained from Schmidt's thermodynamic analysis. 
However, this leads to an optimistic prediction of engine 
performance, since Schmidts’s analysis neglects the 
thermodynamic and fluid frictional losses. Better esti- 
mates of the (K p ] and (C*) matrices may be evaluated 
from a more Inclusive analysis Inherent in the GUMPS 
(Refs. 13 and 14) and HFAST (Ref. 15) design codes, or 


test data. 

The pressure terms P e and P^ are obtainable 
from the design codes, for given piston and displacer 
positions and speeds. For any Stirling engine, the 
partials In EQ (5) are constant with respect to the engine 
dynamics and. thus, can be frozen In analytic computa- 
tions. Additional simplifying assumptions are that the 

partials fjl. £^2 and are all negligible. 

W 0 W r 

Substitution of the reduced EQ (5) into EQ (4) results In 
EQ (7), In 
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which the pressure, position and speed terms are com- 
plex, and can be obtained from the aforementioned 
design codes or test. Equation (7) yields a total of four 
real and imaginary equations from which the real-only 
partials are solved. Subsequently, EQS (8) are obtained 
from EQ (6) and the resulting EQ (5). The expanded 
dynamic EQ (9) results from EQS (1) and (8). 
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Thermoelectric power generation requires connection 
of the FPSE to an electromechanical device, such as a 
linear alternator and an associated load. Modelling of the 
alternator-load system Is discussed next 


ELECTRICAL EQUATIONS OF LA-LOAD MODEL 


An equivalent circuit diagram of the LA is shown in 
Fig. 3. The alternator output terminals are connected 
through a series tuning capacitor C, to a parallel 
combination of a parasitic load R* and an external 
series R l - L*. - C L static load. The capacitor C T serves 
to ensure an electrical resonance within the circuit and, 
hence, a maximum power transfer from the engine- 
alternator system to the connected load. The stability of 
the system may be enhanced by ensuring a coincidence 
or near coincidence of the electrical resonant and the 
mechanical operating frequencies. The other parameters 
In the figure are the stator resistance, R,, the leakage 
Inductance, L,, the magnetizing inductance, L„, and the 
eddy current and hystersis loss term Rg In the magnetic 
core of the LA. 

Application of Kirchhoffs voltage law to loops 1 
through 3 yields EQS (10) to (12) in which the term / 

denotes and. 


(R. ♦ R)i t - RJ l * (L, * LJi, - LJ, ♦ V'c * •, (10) 


RJ\ ~ Lj t * lJ, • 0 <"> 

-R,l, * {R, ♦ ♦ U ♦ - ° < 12 > 

in terms of the engine piston-induced flux change, the 
generated alternator voltage Is: 



(13) 
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The capacitor voltages are stated in EQ5 (14) 
and (15). 


*c * yc r 0*> 

^ - ijc, < 15 > 


Substitution of the state variable EQ (19) into the FPSE 
dynamic equations (9). and into the LA-load equations 
(14) to (18) yields the state-space model of EQ (20), 


• MIX) - (AM 

l VI - iqixi . [qtuj 


(20) 


The term <t> is the flux linking the N turns of the 
LA magnetic material. Solution of EQS (10) to (13) 
results in the derivative currents summarized in EQS (1 6) 
to (18). 
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STATE-SPACE REPRESENTATION 

The external force (EQ (9)) on the piston is 
represented by the term ^v££jl ( . The system state 
variable vector Is defined in EQ (19), 

[XT - [ X'X'XAw v c v c ) - < 19 > 


where the superscript T* denotes transposition of the 
vector. 


The elements of the 9*9 (A]-matrix of EQ (21) are 
defined earlier. The selected state variables are all 
physically measurable and can, therefore, constitute the 
system outputs (Y]. Thus, [Y] is 9*1 vector and the con- 
trol matrix |C) equals a 9*9 Identity matrix (I). If the 
engine oscillation Is initiated by a unity step input to the 
piston, with subsequent motion of the displacer, then the 
input vector (UJ T - (U D U P | « (0 1J. Hence, [D] is a 9*2 
null matrix. 

The matrix [A] contains all the individual and coupled 
variables of the study system. The damping and stiffness 
coefficients contain the engine nonlinear contributions. 
However, the matrices [A], (BJ, [C] and (DJ are an numeri- 
cally constant. Such a linear time-invariant (ITI) system 
can be subjected to small perturbation analysis about an 
operating point, to determine the overall system dynamic 
behavior. 

Evaluation of the eigenvalues of the [A]-matrix 
requires knowledge of all its elements. The only available 
set of complete data pertains to the 5PRE/LA system 
connected to an R-Lload. The static load Is parasitic and 
mainly resistive with small but finite Inductance. There is 
no separate external load as such. Furthermore, the 
existing data shows no core loss resistance for the 
alternator model. Imposition of these simplifications on 
Fig. 3 results in the state variable vector and system 
matrix as defined in EQS (22) and (23), where the 
subscript S denotes simplified electrical model. These 
equations are used in the analysis of the SPRE/LA-load 
system. 
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METHOD OF ANALYSIS 

The SSM of the coupled equations are simulated 
using the MATLAB software. Variations are enforced in 
key system and operating parameters to induce a migra- 
tion of the system eigenvalues in the frequency plane. 
The effect of a parameter on the eigenvalues is deter- 
mined by varying only that parameter about Its nominal or 
design value, while keeping all other parameters at their 
nominal values. A plot of the eigenvalue movement due 
to changes in a given parameter yields the corresponding 
root locus. 

Dynamic analysis derives Its strength from its ability 
to predict the time-domain behavior of an LTI system via 
simple calculations In the frequency domain. This Is 
exemplified by Fig. 4 (Ref. 16) which shows time histories 
for corresponding eigenvalues in the frequency domain. 
The upper part of the plot is typically mirrored in the lower 
half. Exponentially Increasing oscillations are evident for 
the eigenvalue in the right-hand plane. This dynamic 
instability, in the case of a FPSE, will be equivalent to a 
piston overstroke. An eigenvalue on the Jw-axls is 
critically stable, since the oscillation amplitudes are 
constant. This is characteristic of, and desired for a 
FPSE as an oscillator. An eigenvalue in the left-hand 
plane Is dynamically stable. However, a FPSE in this 
mode of operation is described as stalling or falling out, 
since the oscillations decay with time. 

The foregoing discussions are used to characterize 
the simulation results. 

DISCUSSION OF RESULTS 

The interactions between the FPSE/LA subsystem 
and the connected load are Illustrated Figs. 5 to 12. The 
nominal and/or design values of the parameters used in 
the simulation are summarized in Table 1 which repre- 
sents the 'base* or reference system. All the values 


shown are measured quantities. However, zero values 
are enforced for the partials of the heat exchanger 
pressure with respect to the piston and displacer veloci- 
ties. The difficulty associated with measurement of these 
partials casts doubt on the accuracy of their values. 
Unless otherwise indicated, selected parameters are 
varied from 50 to 150 percent of their nominal value, in 
Increments of 25 percent. The arrows connecting the 
*X’s" in the plots signify the directions of the eigenvalue 
movements, due to parametric changes. 

BASE SYSTEM - Figure 5 shows a set of eigen- 
values for the coupled FPSE/LA-load system, and another 
set for the decoupled FPSE and LA-load subsystems. 
The electrical and mechanical roots for the coupled base 
system are identified as E, and MP and MD, respectively. 
The terms MP and MD refer to the piston and displacer 
eigenvalues, respectively. The terms EU, and MPU and 
MDU denote the electrical and mechanical roots for the 
uncoupled FPSE and LA-load subsystems. The uncoupl- 
ing is realized by setting the A (5,4) and A (4,5) elements 
to zero In EQ (23) of the simplified electrical model. 

In the case of the uncoupled subsystems, the com- 
puted 98.6 Hz of the electrical root oompares nearly 
exactly with hand-calculated 98.3 Hz from the characteris- 
tic equation of the electrical subsystem. The average of 
the predicted piston and displacer frequencies is 107.1 
Hz. This is only 2 peroent and 5.8 percent above the 
design and operating frequencies 105 and 101.1 Hz, 
respectively, of the SPRE. It Is noteworthy that the piston 
root MPU in the right-hand plane (RHP) is characteristic 
of an oscillator with exponentially Increasing amplitudes 
towards its mechanical limits. This is consistent with 
Redlich and Berchowltz’s observation (Ref. 2). 

The coupled system shows a slight Increase in the 
average engine frequency, namely, 111 Hz. Again, this 
is within 10 percent of either the design or operating 
frequency of the SPRE. Also, the coupled engine tends 
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to reduce the margin of stability of the electrical system. 
However, the effect of the LA-load subsystem is to push 
the MPU root to the MP position, nearly on the jw-axis. 
Hence, the electrical system acts as an external means 
of control, and forces the engine to perform as an 
oscillator. The engine is said to be critically stable in that 
it oscillates with constant amplitudes. Finally, the unequal 
piston and displacer frequencies Indicate that in the 
practical engine, the piston stored energy/energy dissipat- 
ed in a cycle is not necessarily the same as that of the 
displacer. 

The following discussions will highlight the influence 
of various parameters on the coupled system roots E, MP 
and MD and, hence, the system dynamic stability. 

LA Resistance (R # ) - Varying the alternator resis- 
tance from 50 to 150 percent of its nominal value has 
insignificant effect on the system roots which remain 
almost stationary in Fig. 6. The relatively small value of 
the nominal resistance of 0.0820 may account for its 
ineffectiveness. 

LA Leakage Inductance (LJ - A definite effect on the 
system roots caused by varying L, is evident in Fig. 7. 
The symbols ‘X’ of the base system are identified by self- 
directed arrows from MD, MP and E. Increasing L, from 
50 percent of nominal can place the piston root MP on 
the jw-axis, while improving the damping of the electrical 
root E. The displacer root gets pushed towards the jw- 
axis, but remains in the LHP. Further increase in L, 
beyond the nominal value may destabilize the engine, 
unless this is prevented by some external feedback 
control system. Thus, the nominal L, is considered near 
optimum for the piston root. 

Tuning Capacitance (C T ) - Figure 8 shows that a 50 
to 1 50 percent variation of the nominal tuning capacitance 
has an effect similar to that caused by the leakage induc- 
tance. Here, the optimum capacitance for placing the pis- 
ton root on the Jw-axis Is between 100 and 112 percent 
of the nominal value. 

L oad Resistance (RJ - The pronounced impact of 
the load resistance variation from 50 to 300 percent of its 
nominal value is illustrated in Fig. 9. Increasing Rl 
significantly improves the electrical damping, with only a 
modest change in its frequency. The piston root can not 
only be moved past the jw-axis, but also, its frequency 
decreases gradually. Initially, the frequency of the 
displacer root reduces only slightly, while Its damping 
increases. Beyond 150 percent of the nominal l\, the 
damping reduces with further increase in R t . Examination 
of Fig 9 shows that, as R L increases further, the piston 
root is pulled into the RHP, while the displacer root 
moves towards the jw-axis. These counter-balancing 
effects will tend to stabilize the engine. Since It is the 
larger of the two resistances in the system, the load 
resistance can mask the Influence of the LA resistance. 


The above observation shows that, generally, a 
properly designed parasitic load may be used to force the 
FPSE/LA subsystem to operate as an oscillator; that is, 
one eigenvalue pair on the jw-axis, and all others in the 
LHP. 

Short Circuit Condition - A potentially hazardous 
operating condition in a FPSE/LA-load system is a short 
circuit of the load. This represents an extreme case of 
load impedance variation. Figure 1 0 depicts the positions 
of the system roots for such a condition. The electrical 
root nearly collapses on the jw-axis, without substantial 
change in frequency, when compared to the base system. 
With only the small alternator resistance remaining in the 
electrical circuit, the shorted electrical circuit approximates 
an oscillator. The piston root is forced into the RHP, 
signifying exponentially increasing oscillations. Hence, 
the effect of a sudden short circuit condition is to render 
the engine unstable, with a possible consequence of 
piston overstroke. In the physical system, this situation 
may result from generated energy which is dissipated in 
a minimum Impedance of the electrical system. 

Open Circuit Condition - Another possible but undesir- 
able operating scenario is an open circuit condition. This, 
also, is an extreme form of load inpedance variation. The 
electrical roots in Fig. 11 collapse onto the origin, with 
their damping capability eliminated. The engine becomes 
unstable, as the piston roots are in the RHP. This 
condition Is also possible in a practical system in which 
the generated energy, in the form of power, abruptly has 
no path to flow. 

CONCLUSIONS 

This paper uses the state-space technique to deter- 
mine the effects of parametric variations on the dynamic 
stability of a FPSE/LA-load system. The mathematical 
formulation includes the major thermodynamic effects of 
the engine. The following conclusions are based on the 
evaluation of stability of the SPRE/LA-load system. 
However, the analytical approach and trends of the 
results are applicable to other FPSE/LA-load systems 

The results confirm Redlick and Berchowitz's classical 
control theory-based result that there exists only one 
mode of oscillation for the FPSE. This mode Is denoted 
by a pair of mechanical roots in the right side of the 
complex plane. Such roots must be forced by a control 
mechanism to move onto the jw-axis, to ensure the 
engine behaviour as a constant amplitude oscillator. 

For a given FPSE/LA subsystem, the operating 
frequency can be load dependent. However, the effect of 
the constituent load components depends on the their 
relative magnitudes. A system-designed load generally 
aids in stabilizing the system. Dynamic stability is 
reasonably assured, if the electrical system parameters 
do not change significantly from their design values. For 
the system studied, the parameters with the most impact 
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on the system stability are the LA leakage inductance, the 
tuning capacitance and the load resistance. 

Abnormal operating conditions, such as a sudden 
short or open circuit of the load, can destabilize the 
engine, unless ameliorating controls are in place. 

Generally, the effects of parametric variations on the 
system dynamic stability give an Insight Into the engine- 
altemator-load interactions. This Information can be 
valuable during design stage, and with respect to perfor- 
mance prediction to guide system evaluation in experi- 
mental work. 
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Figure 1 - Block diagram representation of study system. 
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Fig 2 Schematic of a single cylinder 
Free-Piston Stirling Engine 
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Fig. 1 1 Plot of system eigenvalues with open circuit 



Table 1 Summary of System Parameters 


Fig. 9 Plot of system eigenvalues with Rl variation 
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